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Theory of the isotope shift of the centroid energies of light few-electron atoms is reviewed. Nu¬ 
merical results are presented for the isotope shift of the 2 3 P-2 3 S and 2 1 S-2 3 S transition energies 
of 3 He and 4 He. By comparing theoretical predictions for the isotope shift with the experimental 
results, the difference of the squares of the nuclear charge radii of 3 He and 4 He, 8R 2 , is determined 
with high accuracy. 


I. INTRODUCTION 


Helium is the simplest few-electron atom that can be 
calculated ab initio [l| and measured with a very high ac¬ 
curacy. Optical measurements in 3 He 1 and 4 He 0 
has recently reached the relative precision of few parts in 
10 -12 . Such measurements are sensitive to the nuclear 
size effect on the level of 10 -4 . They are approaching the 
level that might allow an improved determination of the 
Rydberg constant, which is currently considered as the 
one of the best known fundamental constants [Bj. The 
accuracy of the present-day experiments is high enough 
for an accurate spectroscopic determination of the nu¬ 
clear charge radii of 3 He and 4 He. Such determination 
would be of particular interest in view of the discrepancy 
for the proton charge radius observed in the muonic hy¬ 
drogen experiment Jf| 0] • The follow-up experiment on 
the muonic helium [8| might soon bring the charge radius 
of helium into focus. 

Unfortunately, the present-day theory of the energy 
levels in helium 0, [HI is not good enough to provide the 
nuclear charge radii of 3 He and 4 He with the required 
precision. A more feasible task at present is the high- 
precision determination of the nuclear charge radii dif¬ 
ference of these isotopes, which is extracted from the iso¬ 
tope shift of the transition energies. In the current paper 
we make a review of the theory of the isotope shift in 
light few-electron atoms, describe our calculation of the 
isotope shift for helium, and determine the difference of 
the squares of the nuclear charge radii of 3 He and 4 He, 
SR 2 . The results of the SR 2 determination have already 
been published in Ref. Q. In this work we present details 
of the calculation and the underlying theory. 

We define the isotope shift of two isotopes as the dif¬ 
ference of the centroid energies. More specifically, we are 
presently interested in the centroid energies of the 2 P 
and 2 S states, which are defined as an average over all 
fine and hyperfine energy sublevels, 


E(2 2S+1 L) 


E. 7 , f ( 2F + 1)U(2 2S + 1 L j , f ) 

(2 / + 1) (2 S' + 1) (2 A + 1) ’ ^ ’ 


where 2S+1 L denotes the state with the angular momen¬ 
tum L and the spin S. This definition of the isotope 


shift differs from the definition sometimes used previously 
(e.g., in Refs, [lj, l2j, |Tl() by the fact that we average out 
not only the hyperfine but also the fine-structure split¬ 
ting. The advantage of using the centroid energy is that 
its theory is much more transparent than that of the 
individual sublevels 2 2S+1 Aj iF . Moreover, the centroid 
transition energy can be directly accessed in an exper¬ 
iment, by measuring all independent fine and hyperfine 
transitions and averaging them out. 

Theory of the isotope shift is much simpler that the¬ 
ory of the energy levels because only very restricted set 
of operators contribute to it. Specifically, only the oper¬ 
ators that depend on the nuclear mass or the nuclear size 
contribute to the isotope shift. In addition, due to the 
averaging over the fine and hyperfine energy sublevels in 
©, all operators that depend the nuclear or the elec¬ 
tron spin do not (almost) contribute to the isotope shift. 
More exactly, such operators contribute only through the 
second-order perturbations, most of which are negligible 
at the level of the present interest. 


II. THEORY 

Within the nonrelativistic QED approach, energy lev¬ 
els of atoms are represented by an expansion in powers 
of the fine-structure constant a , 

E(a) = U (2) + A (4) + U (5) + U (6) + U (7) + ..., (2) 


where E^ = ma n £ l ' n ' > is a contribution of order a n and 
may include powers of In a. Each of is in turn ex¬ 
panded in powers of the electron-to-nucleus mass ratio 
m/M 


( 3 ) 

where denotes the correction of first order in m/M 
and is the second-order correction. For the descrip¬ 
tion of the isotope shift, we need only the part that de¬ 
pends on the nuclear mass M , i.e., corrections £$, £^l, 
etc. 
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A. Nonrelativistic recoil 

The leading term in the expansion (j2Jl is the nonrela¬ 
tivistic energy E^ 2 \ which is the eigenvalue of the non¬ 
relativistic Hamiltonian, 



it up to the third order. The perturbation corrections 
are 



where Z and P are the nuclear charge and momentum, 
respectively; the indexes a and b numerate the electrons, 
p a is the momentum of the electron a, r a is the distance 
between the electron a and the nucleus and r ab is the 
radial distance between the electrons a and b. 

Assuming that the total momentum is zero, 


A 2 ) 

3 — 



_-_ffW 

{So-Ho)' rec 


{So - Ho)' 


ff (2) 



---//vd\ 

{So - H 0 )' 2 rec /J 
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Pa'Pb, 
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a<6 


and rescaling momenta and position operators in units 
of the reduced mass p = mM/(m + M), p a —> pp a and 

/o\ 

f a f a /p, the Hamiltonian H y M ’ is represented as a sum 
of two terms, 


with 


rr(2) _ „ rr(2) , P_ rr(2) 

n M ~ f n oo ' fl rec> 




and 



m 


a<.b 


( 6 ) 

(7) 

( 8 ) 


( 2 ) ( 2 ) 

The eigenvalue of Ho 0 is £/o , the nonrelativistic energy 
in the nonrecoil limit. In the present work we are in¬ 
terested in the recoil (i.e., M-dependent) corrections to 
E^. The correction coming from p is often termed 
as the normal mass shift. It is given by 


where Hq = Ho a is the nonrelativistic Hamiltonian in 

C 2 ) 

the nonrecoil limit and £q = E/o’ is its eigenvalue. We 
checked that, for P states, the perturbation expansion 
results are in a very good agreement with the nonpertu- 
bative treatment. 


B. Relativistic recoil 


The leading relativistic correction to the energy is 
given by the expectation value of the Breit-Pauli Hamil¬ 
tonian H [njj. The Breit-Pauli Hamiltonian with the 
full reduced-mass dependence is given by: 


p p* nZ , , Z _ f a 

+ ~2~ S + ~i Ga ' ^3 X Pa 
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(13) 



( 2 ) 

The second recoil correction is induced by H / ec , it is 
known as the specific mass shift or the mass polariza¬ 
tion. This correction can, in principle, be calculated nu- 
merically to all orders in m/M by including H/ e / into 
the nonrelativistic Hamiltonian, determining the corre¬ 
sponding eigenvalue and subtracting the nonrecoil limit 
M —> oo. This approach works very well for P states 
but becomes numerically unstable for S states (with the 
basis of the fully correlated exponential functions used in 
the present work). Because of this, we obtain the specific 
mass shift corrections by perturbation theory, extending 


where a a is the vector of Pauli a matrices acting on the 
a electron. The Breit Hamiltonian without the reduced 
mass dependence (i.e., with p m) will be denoted by 

4 ) 

' 7 oo • 

We are presently interested in the recoil corrections 
to the Breit contribution. First of all, there is a recoil 
correction induced by the reduced mass in Eq. m, 

S { m]o = (h {4) - Hg>) . (14) 

The first-order recoil correction E^ l is induced by 
the recoil addition to the Breit Hamiltonian and the 
second-order perturbation correction of the non-recoil 
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Breit Hamiltonian H& and the mass-polarization op¬ 


erator H , 


( 2 ) 


4h = (£) («1S) +2 (£)(»£> 


1 


■ ff (2) 


(< £o-H 0 y rec / ’ 

(15) 


C. QED recoil 


( 5 ) 

The leading QED correction £^o in the nonrecoil limit 
is given by [IS 


where the recoil addition to the Breit Hamiltonian is 


(4) = 

rec 9.M 


ab 


“a X p b * ( 7 a ~ P a 


__ ’ CL 1 CL 

r ' r 3 

’a 1 r, 


(16) 


The first term in the above formula represents the 
electron-nucleus spin-orbit interaction, whereas the sec¬ 
ond term is induced by the electron-nucleus orbit-orbit 
interaction. 

The second-order recoil correction £$2 consists of sev¬ 
eral parts, 


c(4) _ c(4) 


’(4) 


7(4) 


M ,2 — ‘-'M,2a ' L 'M,2b ~ r L 'M,2c ' 


(17) 


The first one is the third-order perturbation correction 
of the non-recoil Breit Hamiltonian and two mass- 
polarization operators ITrec, 

c(4) — o/g( 4 ) ^ g( 2 l ^ ff( 2 )\ 

M ’ 2a \ °° (£o-H 0 y iec (£ 0 ~H 0 y rec / 

+ ( H ™ (£ 0 - Ho) 1 (£o - H 0 y H ™) 

- (Hjg) (HQ 

(18) 

The second part is the second-order perturbation correc¬ 
tion of the recoil Breit Hamiltonian Hrtl and the mass- 
polarization operators H Ta [-, 


£ (5) = Ifi 

oo g 


In (Za) 2 + ^ - In ( 


14, , 164 
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where the singular operator r~ 3 is defined by 

= lim J d 3 r <fi* (f) cj>(r) 

0 (r — a) + 4 7 r <5 3 (r) (7 + In a) 


£<* 3 (r«)> 


( 21 ) 


( 22 ) 


where 7 is the Euler constant. The Bethe logarithm is 
defined as 


ln(fc 0 ) = 


J2aPa( H 0-£o) In [2 (Ho — £<])] EfcPb 


2itZ (T.C S3 ( r c) 


(5) 

The recoil correction £ M consists of three parts 
(A + + £ 3 ) j 


(23) 


(24) 


where 

fr = -3^ } + ^S> 3 (r.)>’ ( 25 ) 

a a <6 


£2 = 


2 2 62 8 f ko \ \. r3 , ., 

, 1 -' ' 2^^ J (r a )> 


-_h , ( Za) + _-_i n ^_ 


n. ' u 


(26) 


and ( m/M ) £3 is the first-order perturbation of due 

( 2 ) 

to the mass-polarization operator 7/ r l e c - 


y( 4 ) =2 / _-_77(2)\ (iqi 

c M,2b — z rec ^ _ jj^y CI rec ; . 

The third part is the Darwin term, which reads for the 
1 = 1/2 nuclei, 

^■-(£)’(£ ?'<'•>)■ (20 > 


D. Higher-order recoil corrections 


The QED correction of order m 2 a 6 /M was estimated 
by using the hydrogen results @, 



The Darwin term is zero for spinless nuclei. 
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The first term in the brackets in the above equation is 
due to the reduced mass scaling of the A 50 one-loop QED 
contribution, the second term is the radiative recoil cor¬ 
rection, and the third term is the pure recoil contribution. 

Another important correction of order m 2 a e /M comes 
from the second-order hyperfine correction. We keep only 
the important part of this correction for the 2 P states, 
which is enhanced by the small 2 3 P — 2 1 P energy differ¬ 
ence in the denominator. This correction is referred to 
as the hyperfine mixing contribution and is given by, for 
the 2 3 P state, 


- C/ M,hfsmix 


(2 3 p\H$yp)\ 

Eq(2 3 P) ~ Eq(2 1 P) ’ 


(28) 


where is the operator responsible for the leading- 
order hyperfine structure splitting. In relativistic units, 


H ht! = 9 fJ’N I- 


a *- 


8lT \ - 
— o[r a )(j a 


+ 


3 r a ( a a 


■ r a - (Tq r„ 


X Pa 


(29) 


where ps = e/(2 to) is the Bohr magneton, = 
e/{2m p ) is the nuclear magneton (with m v being the pro¬ 
ton mass), g is the nuclear g factor, and / is the nuclear 
spin. We note that the last term in Eq. (l29l) does not 
contribute to the hyperfine mixing. For the 2 3 P state, 
the hyperfine mixing contribution is given by the same 
expression (1251) with the opposite overall sign. 


E. Nuclear polarizability 


where R is the root-mean-square (rms) nuclear charge 
radius R = \/{R 2 ), A c is the Compton wavelength di¬ 
vided by 2-7 t and ff s is the nonlogarithmic relativistic and 
higher-order correction. The second term in the brackets 
in Eq. m is the leading logarithmic relativistic correc¬ 
tion, known from the hydrogen theory [Bj. It is state- 
independent and can be directly generalized to the he¬ 
lium case. 

The nonlogarithmic relativistic correction ff s is state 
dependent and, moreover, dependent on the nuclear 
model. In the present work, we estimate /f s basing on 
the theory of He + , where it can be obtained by a direct 
numerical evaluation [l8j]. Our calculations, performed 
with quadruple-precision arithmetics for the Is state of 
He + , yield /f s = —0.04 for the Gauss nuclear model, 
/f s = 0.02 for the sphere nuclear model, and /f s = —0.14 
for the exponential-distribution nuclear model. Basing 
on these results, we estimate the corresponding correc¬ 
tion for helium as /f s = 0.0(2). The error ascribed to 
/f s corresponds to the fractional error of 4 x 10~ 5 in Ff s , 
which is negligible at the level of the present experimental 
interest. 

The correction beyond the mean-square charge radius 
Ez comes from the two-photon exchange [h|, which is, 
in relativistic units, 



Since it is proportional to (R 3 ) it depends on the charge 
distribution. However, this effect is almost negligible for 
light atoms, namely the ratio 

Ez Za (R 3 ) _ 

Ef s 2 A C (R 2 ) [ } 


The nuclear polarizability correction to the energy lev- for He amounts to about 4 x 10 -5 . 
els can be represented as 


4oi = -E^ 3 ( r “)> s p°i- ( 3 °) 

a 


III. THEORETICAL RESULTS 


The coefficient 5 po i in the above formula was evalu¬ 
ated [3 as 5 P oi( 4 He) = 2.07 (20) fm 3 and 5 po i( 3 He) = 
3.56(36) fm 3 . Somewhat smaller values were later ob¬ 
tained in Ref. E2 , the difference being not important at 
the level of interest of the present investigation. 


F. Finite nuclear size 


The leading finite nuclear size (fs) correction to the 
centroid energy is, in relativistic units, 


2tt 

Ef s = — m. 


MS’£<£*>.>) 


vm/ A; 
\2 


1 - (ZaY \n(ZaR/X c ) + ( Za) 2 f is 


(31) 


Our main goal in the present investigation is the deter¬ 
mination of the difference of the nuclear charge radii from 
the isotope shift. To this end, we separate our theory of 
the isotope shift into two parts. The first part is the 
theory the isotope shift for the point-like nucleus. The 
second part is the finite nuclear size effect, which is pa¬ 
rameterized as SEf s = C R 2 , where R is the rms nuclear 
charge radius and the coefficient C is calculated numeri¬ 
cally. Now, a comparison of the experimental value with 
the theoretical results for the isotope shift allows us to 
extract the difference of the squared rms nuclear charge 
radii of the two isotopes. 

Numerical results of our calculation of the isotope 
shift of the 2 3 P-2 3 S and 2 1 S'-2 3 S' transitions for the 
point nucleus are presented in Tables Q] and El respec¬ 
tively. The following masses of the helion and 
the alpha particle m a were used in the calculations 
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TABLE I: 3 He- 4 He isotope shift of the 2 1 S — 2 3 5 transition, 
for the point-like nucleus, in kHz. 


pa 2 

-8632 567.86 

MU 2 (m/M) 

608 175.58 

pa 2 (m/M ) 2 

-7319.80 

pa 2 (m/M ) 3 

0.30 

M a 4 

-8954.22 

M a 4 (/V M ) 

6 458.23 

fia 4 (m/M ) 2 

1.84 

m a 5 ( m/M ) 

56.61 

m a 6 (m/M) 

2.75 (69) 

Hfsmix 

80.72 

NPol 

0.20 (2) 

Sum 

-8 034 065.66 (69) 

Ref. [20]“ 

-8034 067.8(1.1) 


“ Corrected by adding the triplet-singlet hfs mixing 
contribution, in order to compensate the difference in the 
definitions of the isotope shift. 


0, m/rrih = 1.819 543 0761 (17) x 1CT 4 and m/m a = 
1.370 933 555 78(55) x 10" 4 . The first four lines of the 
tables display the numerical results for the nonrelativis- 
tic (~ a 2 ) recoil corrections given by Eqs. ©, ®, CD. 
and m, respectively. The next three lines of the ta¬ 
bles contain the results for the relativistic (~ a 4 ) re¬ 
coil corrections given by Eqs. m, (usd, and (HI)-®, 
respectively. The eighth line of the tables contain the 
QED (~ a 5 ) recoil corrections given by Eqs. (El)-®. 
The line labeled as ma 6 m/M contains the hydrogenic 
approximation of the higher-order QED effects as given 
by Eq. ®. The uncertainty ascribed to this correction 
represents our estimation of the uncalculated QED ef¬ 
fects of order a 6 and higher. Finally, the lines labelled as 
“Hfsmix” and “NPol” contain the hyperfine mixing cor¬ 
rection defined by Eq. (1281) and the nuclear polarizability 
correction given by Eq. ®, respectively. 

The uncertainty of the total theoretical values of the 
isotope shift for the point-like nucleus comes from the 
uncalculated higher-order QED effects. Note that our 
estimation of this uncertainty is much larger than the 
one reported previously in Ref. [l]. 

The difference of the theoretical point-nucleus values 
and the experimental isotope shift results should come 
solely from the finite nuclear size effect. The finite nu¬ 
clear size correction to the centroid energy can be param¬ 
eterized as 

Ef s = C R 2 , (34) 

with C being a coefficient to be calculated numerically. 
According Eq. (1311) . the coefficient C has some depen¬ 
dence on R, but it is very weak and can be neglected at 
the level of present interest. 

The results of our calculations of the coefficient C for 
the 3 He- 4 He isotope shift of the 2 3 P-2 3 S and 2 1 S'-2 3 S' 


TABLE II: 3 He- 4 He isotope shift of the 2 3 S — 2 3 P transition, 
for the point-like nucleus, in kHz. 


pa 2 

12 412 458.1 

M a 2 ( fi/M ) 

21243 041.3 

pa 2 (p/M) 2 

13 874.6 

M a 2 (m/M) 3 

4.6 

pa 4 

17 872.8 

m« 4 (m/m) 

-20 082.4 

MO 4 (m/M) 2 

-3.0 

m a 5 (m/M) 

-60.7 

m a 6 (m/M) 

-15.5(3.9) 

Hfsmix 

54.6 

NPol 

-1.1(1) 

Sum 

33 667143.2(3.9) 

Ref. [If 

33 667146.2(7) 


“ Corrected by adding the triplet-singlet hfs mixing 
contribution, in order to compensate the difference in the 
definitions of the isotope shift. 


transition (centroid) energies are 

C{2 3 P-2 3 S) = -1212.2(1) kHz/fm 2 , (35) 
C(2 1 S — 2 3 S) = -214.66(2) kHz/fm 2 . (36) 

These results can be compared with the previously re¬ 
ported values C(2 3 P — 2 3 5') = —1209.8 0 and C(2 l S — 
2 3 S) = -214.40 HU. 


IV. NUCLEAR CHARGE RADII DIFFERENCE 

We now turn to the determination of the nuclear charge 
radii difference from the isotope shift. Specifically, by 
comparing theoretical and experimental values of the cen¬ 
troid energies of different transitions in 3 He and 4 He, we 
extract the difference of the squares of the nuclear charge 
radii, 

SR 2 = i? 2 ( 3 He) — i? 2 ( 4 He). (37) 

Using the experimental results for the 2 1 S — 2 3 S tran¬ 
sition energies in 3 He and 4 He from Ref. Hi and taking 
into account the experimental hyperfine shift of the 2 3 Si 
state, we obtain 5R 2 as described in Table Hill The results 
is 

SR 2 [Rooij 2011] = 1.028 (11) fm 2 . (38) 

This value is somewhat different from the original result 
of Ref. 0, SR 2 = 1.019 (11) fm 2 because of the updated 
theoretical prediction for the 2 1 S' —2 3 S' isotope shift (see 
Table U). 

The second experiment on the helium isotope shift 
suitable for the determination of SR 2 is the measure¬ 
ment of the 2 3 P-2 3 S transition energies in 3 He and 4 He 
in Refs. a a. Using these experimental values and the 
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theory summarized in Table El we obtain (see Table El 
for details) 

(5i? 2 [Cancio 2012] = 1.074 (4) fm 2 , (39) 

which deviates from the result ( 1351 ) by about 4<r. 

The third experiment that can be used for the deter¬ 
mination of SR 2 with comparable accuracy is the mea¬ 
surement of the isotope shift difference F( 3 He, 2 3 P 0 1//2 — 
2 3 £j/ 2 )—P( 4 He, 2 3 P 2 —2 3 Si) performed Ref. [ll|. In this 
experiment, centroid transition energies were not mea¬ 
sured, so further theoretical and experimental input is re¬ 
quired to extract the charge radii difference 5R 2 . In order 
to deduce the experimental value for the isotope shift of 
the 2 3 P-2 3 S centroid energy, we subtract from the exper¬ 
imental result of Ref. [llj] the experimental hyperfme shift 
<5£'hfs(2 3 S'i^ 2 ), theoretical fine shift SE{ s (2 3 P 0 ), and the 
theoretical fine and hyperfme shift 8E\ l t s {2 3 P^ 2 ). Hav¬ 
ing obtained the isotope shift of the 2 3 P-2 3 S' centroid 
energy, we deduce the charge radii difference SR 2 with 
the help of the theory in Table El The whole determina¬ 
tion is summarized in Table [V] The result is 

SR 2 [Shiner 95] = 1.066 (4) fm 2 . (40) 

We observe that the two results from the 2 3 P — 2 3 S 
transitions, HMD and (SOI), are in reasonable (although, 
marginal) agreement with each other, but deviate sig¬ 
nificantly from the result deduced from the 2 1 11 S — 2 3 S 
transition ( 1551 ) . The reason of this discrepancy is not 
known at present. 


V. CONCLUSION 


In the present work we have summarized the theory 
of the isotope shift of the centroid energies of light few- 
electron atoms. Numerical results are reported for the 
isotope shift of the 2 3 P-2 3 S and 2 1 S-2 3 S transition en¬ 
ergies in 3 He and 4 He. As compared to the previous 
evaluations 0,13, the higher-order recoil [^a 2 {n/M) 3 
and ga 4 (/i/M) 2 ] and the nuclear polarizability correc¬ 
tions were accounted for and the higher-order QED ef¬ 
fects [ma 6 (m/M)\ were estimated more carefully. By 
comparing the theoretical predictions with the experi¬ 
mental results for the isotope shift, we have determined 
the the difference of the squares of the nuclear charge 
radii of 3 He and 4 He, SR 2 . The results for SR 2 derived 
from the 2 3 P-2 3 S and 2 1 S-2 3 S transitions are shown to 
disagree with each other by about 4 standard deviations. 
The reason of this disagreement is presently not known. 
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TABLE III: Determination of the charge radii difference SR 2 from the measurement by Rooij et al. in Ref. 2], in kHz. The 
remainder SE is proportional to SR 2 , 5E = C SR 2 , with the coefficient C given by Eq. (1351) , see text. 


P( 3 He, 2 1 s F=1/2 - 2 3 S f=3/2 ) - P( 4 He, 2*5 - 2 3 5) 

-5 787 719.2(2.4) 

Exp. [2] 

SE hts (2 3 S 3 ^ 2 ) 

-2 246 567.059(5) 

Exp. [22, 231 

—SEi so (2 1 S — 2 3 S) (point nucleus) 

8 034 065.66 (69) 

Theory, Table [J 

SE 

-220.6(2.4) 



TABLE IV: Determination of the charge radii difference SR 2 from the measurements by Cancio Pastor et al. in Refs. a a , in 
kHz. The remainder SE is proportional to SR 2 , SE = C SR 2 , with the coefficient C given by Eq. (13611 . see text. 


£( 3 He, 2 3 P - 2S) (centroid) 

276 702 827 204.8 (2.4) 

Exp. [3J 

-£( 4 He, 2 3 P - 25) (centroid) 

-276 736 495 649.5(2.1) 

Exp. [4, 24f 

— SE iBO (2 3 P — 2 3 5) (point nucleus) 

33 667143.2 (3.9) 

Theory, Table IIII 

SE 

-1301.5(5.0) 



“ obtained by combining the 2 3 5 i — 2 3 Po,i transition energies from Ref. [4| and the 2 3 Po — 2 3 P 2 energy from Ref. nr 


TABLE V: Determination of the charge radii difference SR 2 from the measurement by Shiner et al. in 0 , in kHz. The 
remainder SE is proportional to SR 2 , SE = C SR 2 , with the coefficient C given by Eq. ([361) , see text. 


£( 3 He, 2 3 P 0 1/2 - 2 3 S 3 1 ' 2 ) - E( 4 He, 2 3 P 2 - 2 3 5i) 

810 599.0(3.0) 

Exp. [11] 

SE h{s (2 3 S 3 1 /2 ) 

-2 246 567.059(5) 

Exp. [22, 23] 

SE !b (2 3 P 2 ) 

-4 309 074.2(1.7) 

Theory [25] 

—<5£f s , h f s (2 3 P 0 1/2 ) 

-27923 393.7(1.7) 

Theory [26] 

—SEi so (2 3 P — 2 3 S) (point nucleus) 

33 667143.2(3.9) 

Theory, Table IIII 

SE 

-1292.8(5.2) 















